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Abstract: Llight-element abundances probing big bang nucleosynthesis (BBN) and precision data from cosmology 
probing the cosmic microwave background (CMB) decoupling epoch have hinted at the presence of extra 
relativistic degrees of freedom. This is widely referred to as "dark radiation", suggesting the need for new light 
states in the UV completion of the standard model (SM). We provide a brief and concise overview of the current 
observational status of such dark radiation and investigate the interplay between two possible interpretations of the 
extra light states: the right-handed partners of three Dirac neutrinos (which interact with all fermions through the 
exchange of a new heavy vector meson Z') and dark matter (DM) particles that were produced through a non- 
thermal mechanism, such us late time decays of massive relics. This model ties together cosmological indications 
of extra light states and the production of the heavy vector particle at the Large Hadron Collider (LHC). 
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To accommodate new physics in the form of extra rela- 
tivistic degrees of freedom (r.d.o.f.) it is convenient to ac- 
count for the extra contribution to the SM energy densi- 
ty, by normalizing it to that of an "equivalent" neutrino 
species. The number of "equivalent" light neutrino species, 
iV e ff = (pr — py) /pv L < quantifies the total "dark" relativistic 
energy density (including the three left-handed SM neu- 
trinos) in units of the density of a single Weyl neutrino: 
p VL = (77T 2 /120)(4/ll) 4 / 3 :r r 4 , where p y is the energy den- 
sity of photons (with temperature Ty) and Pr is the total 
energy density in relativistic particles fl~|. A selection of the 
most recent measurements of Ng& together with the 1 a con- 
fidence intervals from various combinations of models and 
data sets are shown in Fig. [T] The data hint at the presence 
of an excess AN above SM expectation of A^ff = 3.046 |2|. 

One of the most striking results of the Planck space- 
craft is that the best-fit Hubble constant has the value 
h = 0.674 ±0.012 [310This result deviates by more than 
2a from the value obtained with the Hubble Space Tele- 
scope, h = 0.738 ± 0.024 |4 |. The impact of the new h de- 
termination is particularly complex in the investigation of 
A?eff. Combining CMB observations with data from bary- 
on acoustic oscillations (BAO) the Planck Collaboration 
reported A^ff = 3.30 ±0.27. Adding the Hq measurement 
to the CMB data gives Ngg — 3.62 ± 0.25 and relieves the 
tension between the CMB data and Hq at the expense of 
new neutrino-like physics (at around the 2.3a level). We 
have nothing further to add to this discussion, but it seems 
worthwhile to analyze any other independent set of data |5 1. 

Several models have been proposed to explain AN. To 
accommodate recent Planck observations in this brief com- 
munication we re-examine two of them: models based on 
milli-weak interactions of right-handed partners of three 
Dirac neutrinos |6 1, and models in which the extra relativis- 
tic degrees of freedom are related to possible dark matter 
(DM) candidates produced via decay of heavy relics [7 1 j^] 
An interesting consequence of such a non-thermal DM sce- 
nario is that if the lifetime of the decaying particle is longer 
than about 10 3 seconds, the expansion history of the uni- 
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Fig. 1: A selection of the most recent cosmological N e ff 
measurements and the 1 a confidence intervals from various 
combinations of models and data sets. BBN findings are 
shown in blue and those from the CMB epoch in red. 



verse during the era of BBN will not have DM contributions 
to N e f[. Moreover, if there is a light DM partcile that annihi- 
lates to photons after the Vl have decoupled, the photons 
are heated beyond their usual heating from e annihilation, 
reducing the late time ratio of neutrino and photon tempera- 
tures (and number densities), leading to A^ff < 3 [8]. This 
opens the window for addition of one or more Vr neutrino 
flavors and still be consistent with N e s = 3. 
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We adopt the usual convention of writing the Hubble constant 
at the present day as Hq = 100 h km s~ 1 Mpc~ 1 . For / = today, 
the various energy densities are expressed in units of the critical 
density p c ; e.g., the matter density O^m = pjy[(today)/p c . 
Vs's have long been suspected to contribute to N e g |9|. 
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We begin by first establishing the contribution of right- 
handed neutrinos to A^ff, that is AAV as a function of the 
Vr decoupling temperature. A critical input for AAV is the 
relation between the r.d.o.f. and the temperature of the 
primordial plasma. This relation is complicated because the 
temperature which is of interest for right-handed neutrino 
decoupling from the heat bath may lie in the vicinity of 
the quark-hadron cross-over transition. To connect the 
temperature to an effective number of r.d.o.f. we make use 
of some high statistics lattice simulations of a QCD plasma 
in the hot phase, especially the behavior of the entropy 
during the changeover iflOl . 

Taking into account the isentropic heating of the rest of 
the plasma between 7^ ec and 7^ ec decoupling temperatures 

we obtain AAV = 3[g s (T v d ^)/g s (T^)} 4 / 3 , where g s (T) is 
the effective number of interacting (thermally coupled) 
r.d.o.f. at temperature T; e.g., gs(Ty^ c ) — 43/4. At energies 
above the deconfinement transition towards the quark gluon 
plasma, quarks and gluons are the relevant fields for the 
QCD sector, such that the total number of SM r.d.o.f. is 
g s = 61.75. As the universe cools down, the SM plasma 
transitions to a regime where mesons and baryons are the 
pertinent degrees of freedom. Precisely, the relevant hadrons 
present in this energy regime are pions and charged kaons, 
such that g s — 19.25 ifTTl . This significant reduction in 
the degrees of freedom results from the rapid annihilation 
or decay of any more massive hadrons which may have 
formed during the transition. The quark-hadron crossover 
transition therefore corresponds to a large redistribution of 
entropy into the remaining degrees of freedom. Concretely, 
the effective number of interacting r.d.o.f. in the plasma at 
temperature T is given by g s (T) ~ r{T) (Nr + ~^Np) , with 
Nb =2 for each real vector field and Np = 2 for each spin- 
5 Weyl field. The coefficient r{T) is unity for leptons, two 
for photon contributions, and is the ratio s(T)/s^ for the 
quark-gluon plasma. Here, s(T) (ssb) 1S the actual (ideal 
Stefan-Bolzmann) entropy shown in Fig [2] The entropy rise 
during the confinement-deconfinement changeover can be 
parametrized, for 150 MeV < T < 500 MeV, by 
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For the same energy range, we obtain 

g s (r)~47.5 r(T) + 19.25. 



(1) 



(2) 



In Fig.|2]we show g s (T) as given by (j2j). Our parametriza- 
tion is in very good agreement with the phenomenological 
estimate of lfT2ll . 

If relativistic particles are present that have decoupled 
from the photons, it is necessary to distinguish between 
two kinds of r.d.o.f.: those associated with the total energy 
density g p , and those associated with the total entropy 
density g s . Since the quark-gluon energy density in the 
plasma has a similar T dependence to that of the entropy 
(see Fig. 7 in £[0)), we take g p (T) ~ r(T) (N B + ^ F ) . 

The right-handed neutrino decouples from the plasma 
when its mean free path becomes greater than the Hubble 
radius at that time 
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Fig. 2: Left: The parametrization of the entropy density 
given in Eq. ^ (dashed line) superposed on the result 
from high statistics lattice simulations ifTOll (solid line). 
Right: Comparison of g s (T) obtained using Eq. (|2| (dashed 
line) and the phenomenological estimate of lfl2l (solid line). 
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Fig. 3: The (green) shaded area shows the 1(7 region al- 
lowed from decoupling requirements to accommodate 
Planck + BAO, whereas the (yellow) cross-hatched area 
shows the region allowed from decoupling requirements to 
accommodate Planck + Ho. We have taken ,3tf = 0.5 (left) 
and X = 2.5 (right). 
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chiral states, gj are the chiral couplings of the Z' for the 
6 relevant species (see below), and the cosntant J?T = 
0.5 (2.5) for annihilation (annihilation + scattering) |6|. In 
the second line of (j5ji we set g s ~ g p . 

The physics of interest takes place in the quark gluon 
plasma itself so that we will restrict ourselves to the fol- 
lowing fermionic fields in the visible sector, [3ur] + [3/1r\ + 
[3s*] + [iv L + e L + pL L ] + [e R +li R ] + [3u L + 3d L + 3s L ] + 
[3Vr], and their contribution to g p . This amounts to 28 Weyl 

fields, translating to 56 fermionic r.d.o.f.; Y.f=i = 28. To 
illustrate we calculate g for two candidate models. The first 
is a set of variations on D-brane constructions which do not 
have coupling constant unification, 0.29 < g < 0.39 03). 
The second are U ( 1 ) models which are embedded in a grand 
unified exceptional group, 0.09 < g < 0.46 [9|. 
Substituting Q and Q into Q we obtain 
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In Fig. [3] we show the region of the parameter space 
allowed from decoupling requirements to accommodate 
contributions of AA/ V within the 1 o~ region of Planck data. 
It is important to stress that the LHC experimental limits on 
M z i for null signals for enhancements in dilepton searches 
entail M z > > 2.3 TeV at the 95% CL El, whereas limits 
from dijet final states entail M z > > 4 TeV at 95% CL |fl31 . 

Next, in line with our stated plan, we derive quantitative 
bounds on the fraction of non-thermally produced DM 
particles, which are relativistic at the CMB epoch. The 
assumption herein is that of the total DM density around 
today, Qdm, a small fraction, / = Q,x/Qdm, is of particles 
of type X, produced via decay of a heavy relic X' with mass 
M' and lifetime t: X' -> X + y^At any time after the decay 
of X' the total DM energy density is found to be 
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where M is the mass of the X particle, nx(t) its number 
density, and a{t) is the expansion scale factor normalized by 
a(today) = 1. The scale factor dependence on the Lorentz 
boost is given by 
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where 
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At the present day the X particles are non-relativistic, 
i.e. Clx (today) =M nx (today) /p c . To obtain such a non- 
relativistic limit we demand the magnitude of the second 
term in the expansion of ([8]) to be greater than the third term, 
which results in [a{%) / 'a[t)} 2 [y 2 (t) — 1] < 4. Contrariwise, 
by this criteria the particle X is relativistic if 7(0 > V5- 

To obtain the X component of N e g we assume that X 
particles decouple from the plasma prior to Vl decoupling, 
conserving the ratio Ty/T Vl from SM cosmology, 
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where Py{t E Q) = PcQy/a 4 ( ? eq)- The radiation-matter equal- 
ity time is defined by the condition Pr = Pm, which implies 
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with AW = AN V + ANx- By expressing the scale factor in a 
piece-wise form, 
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Fig. 4: Left: Contours of constant ANx in the R— t/?eq 
plane, with / = 0.01, for the case in which AN V = 0. The 
(green) solid line indicates the upper limit on A^V from 
Planck + BAO. The band between the dashed (yellow) lines 
corresponds to the allowed A^V from Planck + Hq. The 
crosshatched area pertains to the region of the parameter s- 
pace for which the X particles are non-relativistic at the CM- 
B epoch, and should be therefore not taken under consid- 
eration. Right: The shaded regions show BBN constraints 
on the decay X' — » X + y \ 17 1. The various contributions to 
ANx fall on straight lines. 
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Note that the combination ht E Q is independent of h. 

At this point it is worth exploring the quanity 7(?eq) and 
its dependence on h, 

YM - ^ + {a(T)/a(t EQ )} 2 (fl2-l) 2 /(4/? 2 ) 

l + (T/f EQ ) (fl2-l) 2 /(4fl 2 ), (12) 

where R—M'/M and we have assumed that X' decays be- 
fore the time of matter radiation equality. Interestingly, if 
time is scaled to t E Q the analysis is h independent. This is 
important because of the recent tension between astrophysi- 
cal data and cosmological observations. The systematics on 
h will seriously influence the dependence on A^ff- Putting 
all this together ( 10 1 can be rewritten as 
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Note that ( [73) scales with the ratio Q. e ,m/^m — 0.84, which 
does not depend on the Hubble parameter. Finally, solving 
(fT3) for AN X gives 



AN X = (7.4 + AA/ v ) 
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In Fig. [4] we show contours of constant ANx in the R — 
l/t E Q plane, for the case in which AN V = 0. As expected, 

3. We do not at present consider the more complicated scenario in 
which high energy neutrinos are among the decay products 1 16 1. 
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to produce a given ANx contribution, the required ratio of 
masses diminishes with increasing lifetime. 

We now verify that X' — » X + y does not drastically alter 
any of the light elemental abundances synthesized during 
BBN. Following IlTTII . we assume the photons injected 
into the plasma rapidly redistribute their energy through 
scattering off background photons and through inverse 
Compton scattering. As a consequence, the constraints 
from BBN are (almost) independent of the initial energy 
distribution of the injected photons, and are only sensitive 
to the total energy released in the decay process. In the spirit 
of H8 |, we conveniently write the electromagnetic energy 
release as £ y = E y Y X ', where E y = (M' 2 -M 2 )/(2M') is the 
initial electromagnetic energy release in each X' decay and 
Y x i = rix' I n y G is the number density of X' before the decay, 
normalized to the number density of background photons 
riy G — 2 £(2) Ty / n 2 . For Y x i, each X' decay produces one 
X, and so the X' abundance may be expressed in terms of 
the present X abundance through 



} Z'( T ) = Yx,t = Ijf, today = 
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2 . 26xl0 -i4^X^!^ ; (16) 
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GeV. 



The thorough analysis of electromagnetic cascades reported 
in ifTTll reveals that the shaded regions of Fig. |4] are ruled 
out by considerations of light elemental abundances pro- 
duced during BBN. The various regions are disfavored by 
the following conservative criteria: (i) D/H < 10~ 49 (low); 
(ii) D/H > 10~ 4 - 3 (high); (in) 7 Li/H < lO" 10 05 ; (iv) primor- 
dial 4 He abundance < 0.227. The straight lines represent 
several combinations of R and t/?eq producing the ANx 
indicated in the labels. All straight lines intersect the BB- 
N bounds at about log 10 (T//f£g) = —8.2. The constraints 
from BBN are weak for early decays because at early times 
the universe is hot and thus the X' secondary photon spec- 
trum is rapidly thermalized, leaving just a few high-energy 
photons that cannot alter the light elemental abundances. 
However, for T/Yeq > 10~ 8 2 BBN excludes most of the 
relevant parameter space. 

Finally, we combine (|6]) and \\5) to account for both 
AN V and AN X components. In Fig.|5]we show contours of 
constant AN in the M z > /(Kg) vs. (i/t m )(R 2 - l) 2 /(2R) 2 
plane. Interestingly, despite the very restrict limit from 
Planck + B AO data, there is room within the 1 o~ allowed 
iV e ff region to accommodate contributions from both non- 
thermal DM and right-handed neutrinos which interact with 
all fermions via a Z' gauge boson that is within the LHC 
discovery reach. Note that for (say) ANx ~ 0.3, the par- 
ticular range of DM parameters includes (t/?eq ~ 10~ 8 ' 2 , 
X«6x 10 4 ), which is not excluded by BBN considera- 
tions. 
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